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Abstract 

In the second part of investigation of metric nonlinear gravity theories we study 
a fundamental criterion of viability of any gravity theory: the existence of a stable 
ground state solution being either Minkowski, de Sitter or anti-de Sitter space. 
Stability of the ground state is independent of which frame is physical. In general 
a given theory has multiple ground states and splits into independent physical 
sectors. The fact that all L = f{gaf3,Rfiu) gravity theories (except some singular 
cases) are dynamically equivalent to Einstein gravity plus a massive spin-two and 
a massive scalar field allows to investigate the stability problem using methods 
developed in general relativity. These methods can be directly applied to L = f{R) 
theories wherein the spin-two field is absent. Furthermore for these theories which 
have anti-de Sitter space as the ground state we prove a positive-energy theorem 
allowing to define the notion of conserved total gravitational energy in Jordan 
frame (i.e., for the fourth-order equations of motion). As is shown in thirteen 
examples of specific Lagrangians the stability criterion works effectively without 
long computations whenever the curvature of the ground state is determined. An 
infinite number of gravity theories have a stable ground state and further viability 
criteria are necessary. 

PACS numbers: 04.50. +h, 98.80.Jk 

1 Introduction 

In our previous work [1] (hereafter cited as Paper I) we have argued that the cosmolog- 
ical observations based on the spatially flat Robertson-Walker spacetime (usually fitted 
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by the ACDM model) are unsuitable and insufficient to reconstruct the Lagrangian of 
the true gravity theory which correctly accounts for the present epoch of cosmic accel- 
eration. The very fact that the observational data may be fitted by a huge collection 
of diverse Lagrangians clearly indicates that the idea of reconstructing the correct the- 
ory from cosmology is implausible and a deeper investigation con&ms that the task 
is impossible. If one believes that general relativity needs some modifications, these 
should not be directly induced from the approximate data. Rather, as any other new 
physical theory, a new gravity theory should be based on new concepts and ideas. In 
the present case the modifications consist in replacing the Einstein Lagrangian L = R 
by a nonlinear function of Riemann tensor and it is difficult to give a physical idea 
which would uniquely choose the correct function. One should rather apply various 
criteria taken from classical field theory to maximally narrow down the class of viable 
theories. In our opinion the most fundamental condition is the existence of a maximally 
symmetric stable ground state solution to the equations of motion in pure gravity. Or- 
dinary matter cannot destroy stability of a ground state if it is stable against purely 
gravitational perturbations. Stability of the ground state is independent of which frame 
(the set of dynamical field variables) is regarded as physical. Hence the stability can be 
investigated in Einstein frame where the most reliable methods of checking it have been 
developed. 

In the case of L = f{R) the stability is verified in this frame in a very effective and quick 
way. The main step is to solve an algebraic equation to determine the curvature R of 
a ground state. If the equation is solvable for a given Lagrangian one easily establishes 
whether the state is stable. It is worth noting that our method based on the dominant 
energy condition for the scalar component of gravity, applies to generic (inhomogeneous) 
perturbations of a ground state, while most authors assume homogeneous or at most 
spherically symmetric perturbations. 

We emphasize that stability of the ground state and other possible viability condi- 
tions concern physical viability of a gravity theory from the viewpoint of field theory, 
i.e. they concern its internal structure and its relationships to other physical theories. 
There are many viable gravity theories. At the present level of knowledge there is no 
system of selection rules (i.e. viability conditions) allowing one to uniquely determine 
the correct theory. Hence physical viability has a restricted meaning and does not mean 
that a viable theory necessarily ffis in a satisfactory way some empirical data. Actually 
most of viable gravity theories is in conflict with observations. In particular the physical 
viability (which in the present work coincides with the stability of the ground state) is 
independent of a cosmological viability introduced by Amendola et al. [2]. According to 
these authors the cosmological viability means a satisfactory evolution of the universe 
in the fiat {k = 0) Robertson-Walker spacetime: long matter era with the cosmic scale 
factor a oc 

^2/3 

prior to a late-time acceleration epoch. 
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The present work deals with various aspects of the stabihty problem: determination 
of ground state solutions and their multiplicity, notion of stability, total energy and its 
relationship to stability of a ground state, reliability of various methods of checking the 
stability. Finally we formulate a stability condition in terms of a potential in Einstein 
frame for the scalar component of a L = f{R) gravity. We apply the condition to 
thirteen specific Lagrangians, mainly taken from the literature and sho\\0 that some 
theories which are cosmologically viable are physically untenable. 



2 Stability of a maximally symmetric ground state 

A minimal requirement that may be imposed on a gravity theory to be viable is that 
it has a classically stable maximally symmetric ground state solution. In some classical 
field theories, e.g. in Liouville field theory [4] a ground state may not exist, but in grav- 
itational physics the existence of a ground state hardly needs justification. In a metric 
gravity theory gravitational interactions are manifested by the dynamical curvature of 
the spacetime, hence in the absence of these interactions the spacetime should be either 
flat or maximally symmetric with the nongeometric components of the gravitational 
multiplet equal to zero or covariantly constant. Therefore the spacetime of the ground 
state for any NLG theory may be Minkowski, de Sitter (dS) or anti-de Sittecl (AdS) 
space. For simplicity we assume spacetime dimensionality d = 4 although our argu- 
ments (with slight modifications) will also hold in d > 4. Classical stability means that 
the ground state solution is stable against small excitations of the (multicomponent) 
gravitational field and small excitations of a given kind of matter sources, i.e. there are 
no growing in time perturbation modes. In principle a viable classical field theory may 
admit a semiclassical instability: the ground state is separated by a finite barrier from a 
more stable (in the sense of lower energy) state and can decay into it by a semiclassical 
barrier penetration [5j. We shall not consider this possibility and focus our attention on 
classical stability, hereafter named stability. 

A question that may arise at the very beginning of investigation of the problem is 
whether a metric NLG theory, being a higher derivative one, can at all be stable [6]. In 
point particle mechanics one may invoke to this end the old famous Ostrogradski the- 
orem to the effect that if a mechanical Lagrangian depends on second and higher time 
derivatives of the particle positions (which cannot be eliminated by partial integration) 
the corresponding Hamiltonian is linear in at least one canonical momentum and thus is 
unbounded from below. As a consequence there are both positive and negative energy 

^The problem of the cosmic acceleration applying f{R) gravity was also investigated within the 
Kaluza-Klein framework [3] and some results achieved there are akin to ours. 

^By anti-de Sitter space we always mean the covering anti-de Sitter space without closed timclike 
curves. 



3 



states and if the particles are interacting the theory is unstable since any solution de- 
cays explosively due to self-excitation: unlimited amount of energy is transferred from 
negative energy particles to positive energy ones. By analogy, the same (or rather more 
drastic and violent) instability is expected to occur in classical (and quantum) field 
theory with higher time derivatives. Thus a generic NLG theory should be inherently 
unstable and hence unphysical. We admit that the problem is important and deserves a 
detailed investigation. Here we wish only to make short comment on how it is possible 
to avoid this conclusion. 

We stress that the Ostrogradski theorem is a rigorous "no-go theorem" in classical and 
quantum point particle mechanics [6] while in metric NLG theories it may only be con- 
jectured by analogy. In fact, a mechanical Hamiltonian determines energy and if it is 
indefinite (and unbounded from below) it signals that self-excitation processes are likely 
to occur. Recall that a metric theory of gravity is based on the equivalence principle 
what implies that the notion of gravitational energy density makes no sense. Yet in 
a field theory in Minkowski space the field energy density is equal to the Hamiltonian 
density and the latter is (for known fields) positive definite. In the canonical ADM 
formalism in general relativity the canonical momenta are defined in an intricate way 
(including constraints), not akin to that in point mechanics and the total ADM energy 
is to large extent independent of the detailed form of the Hamiltonian density (which 
is indefinite). Therefore in general relativity the relationship between stability (under- 
stood as the positivity of energy, see below) and the form of the Hamiltonian density is 
very indirect, practically broken. In metric NLG theories the Legendre transformations 
from Jordan frame to Helmholtz- Jordan frame (HJF) and Einstein frame (EF) map the 
higher derivative theory to Einstein gravity plus nongeometric components of the multi- 
plet which dynamically act as some matter fields, therefore the stability problem in these 
theories is reduced to that in the latter theory. The Ostrogradski theorem may rather 
serve as a warning that some troubles may appear there and in fact troubles were found 
(the ghost-like behaviour of the massive spin-2 component of gravity) without resorting 
to it. Note that the notion of "inherently unstable theory" is imprecise: stability always 
concerns a given solution. And what is really required from a viable gravity theory is ex- 
istence of a stable ground state solution; stability of excited states is a different problem. 

In the physical literature there is some confusion concerning stability since there are 
actually two notions of stability: dynamical stability (stability of evolution) meaning 
that there are no growing modes and stability as a consequence of positivity of total 
energy. It has been believed for a long time that the two notions are identical and 
since investigations of energy are relatively easier the research was first centered on it. 
Stability in the context of energy was developed in a series of papers which will be 
referred here to as "classical works". Positivity of the ADM energy implies stability of 
Minkowski space. The notion of this energy was then extended to the Abbott-Deser 
(AD) energy for the spacetimes which are asymptotically de Sitter or anti-de Sitter [7]. 
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Applying this notion it was shown that vacuum dS is hnearly stable [7j while AdS space 
is both linearly and nonlinearly stable in vacuum [7j and in presence of any matter sat- 
isfying the dominant energy condition (DEC) in any dimension d > 4 [HI [9] . However it 
was found that stability does not necessarily result from the positivity of energy: there 
are situations in which the positive energy theorem holds and instabilities develop [10]. 
Thus dynamical stability (no growing in time perturbation modes) and positivity of en- 
ergy are quite different unrelated things. Stabihty of evolution requires mathematically 
rigorous investigations. 

In the rigorous approach it was shown that Minkowski space is globally dynamically 
stable: in vacuum [JJJ, in presence of the electromagnetic field [12] or of the linear 
massless scalar field [13] . Vacuum de Sitter space is globally stable in four [H] and any 
larger even number of dimensions [H]. Inclusion of matter is difficult: global stability 
of dS space was proved only in the case of Yang-Mills fields (in d = 4) [16] and for a 
scalar field with a very specific potential [T7\; its stability for all other forms of matter 
is unknown. Even less is rigorously known about stability of anti-de Sitter space: it 
is globally linearization stable for the Maxwell and linear scalar field [TH] and for the 
vacuum case Friedrich [T^ proved finite time nonlinear stability. There are no rigor- 
ous global results, it is only believed that vacuum AdS space is dynamically stable and 
nothing has been investigated in the case of self-interacting scalar fields. 

While the fully reliable rigorous results are quite modest from the standpoint of a 
physicist dealing with gravitational fields generated by a rich variety of matter sources, 
the classical theorems based on the positivity of energy are, from the viewpoint of mathe- 
maticians, of rather little reliability pO]. In proving the rigorous theorems only the exact 
field equations are relevant and the dominant energy condition does not play explicitly 
any role. However in the few cases where matter sources are present, DEC does hold. 
It is therefore reasonable to conjecture that Minkowski, de Sitter and anti-de Sitter 
spaces are globally nonlinearly stable only if any self-gravitating matter does satisfy the 
condition. The conjecture is supported by outcomes found in the linear approximation 
to semiclassical general relativity where the expectation value < 0|T^j/|0 > cannot sat- 
isfy DEC due to the particle creation by the gravitational field. In the presence of the 
electromagnetic, neutrino and massless scalar fields Minkowski space is linearly unstable 
[21] and similarly a minimally coupled quantum scalar field renders Sitter space linearly 
unstable [22]. 

All the aforementioned papers deal with solutions to Einstein field equations. Re- 
cently Faraoni [5^ studied stability of vacuum dS space in restricted NLG theories in 
Jordan frame for the fourth-order field equations. The dS metric can be presented in 
the form of the spatially fiat Robertson- Walker spacetime and he has applied the gauge 
invariant formalism of Bardeen-Ellis-Bruni-Hwang for perturbations of Friedmann cos- 
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mology. The formalism works for any field equations in this background and he proves 
linearization stability of dS space: scalar and tensor metric perturbations are fading or 
oscillating at late times provided the Lagrangian L = f{R) satisfies some inequality. 
In this formalism the physical meaning of this crucial inequality is unclear. It turns 
out that the condition is equivalent to the condition that the (positive) potential for 
the scalar component of gravity in Einstein frame attains minimum at dS space being a 
ground state solution, see sect. 6. The BEBH formalism does not apply to perturbations 
of AdS space since its metric cannot be expressed as the spatially flat R-W spacetime. 
It is interesting to see that in most papers on NLG theories it is assumed that a curved 
ground state is necessarily dS space while AdS space is omitted without mentiorif]. 

We shall investigate stability of the maximally sjTumetric ground state solutions 
in various NLG theories in a coordinate independent manner. We presume that the 
classical works provide the correct assumptions under which the dynamical stability of 
these solutions will be rigorously proved in future. We shall work in Einstein frame where 
the only source for the metric is the scalar field component of gravity since on physical 
grounds it is stability of pure gravity that is crucial. Moreover we argue in section 4 
that inclusion of matter (e.g. perfect fluid) does not affect stability of the solutions. We 
emphasize that stability of a candidate ground state solution is independent of which 
frame is regarded as physical since boundedness of solutions remains unaltered under 
Legendre transformations. The method based on positivity of total ADM or AD energy 
works directly only in Einstein frame. The energy-momentum tensor of the scalar 
satisfies the dominant energy condition if and only if its potential is nonnegative. Thus 
satisfying DEC for the field becomes an effective viability criterion for restricted NLG 
theories. 



3 Candidate ground state solutions 



We shall now investigate existence of candidate ground state (CGS) solutions, i.e., max- 
imally symmetric (dS, AdS or Minkowski space) solutions in a restricted NLG theory 
with L = f{R) for arbitrary /. A CGS solution becomes a true physical ground state 
solution (named vacuum) if it is stable. We assume that the Lagrangian has the same 
dimension as the curvature scalar, [f{R)] = [R] = (length)"^, and the signature is 
( — h ++). The field equations in Jordan frame are 



or 



E,,{g) ^ f'{R)R,. - \f{R)g,u - V,VJ\R) + g,.nf\R) 



f{R)R,u - V,VJ\R) + \g,. \^f{R) - Rf\R) 







(1) 



^AdS space is mentioned as a possible ground state e.g. in 
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here /' = ^ and □ = ^''"^V^Vjy. In general f{R) cannot be everywhere smooth and the 
nonhnear equations (1) require / be piecewise of clas^. 

A CGS solution exists if and only if the field equations (1) admit a class of Einstein 
spaces, R^u = \\g^u, for some curvature scalar i? = A, as solutions. Since A = const 
and assuming that A lies in the interval where /(A), /'(A), /"(A) and /'"(A) are finite, 
eqs. (1) reduce to an algebraic equation 

A/'(A)-2/(A) = 0. (2) 

This equation was first found by Barrow and Ottewill p6] and then rediscovered many 
times. In general this equation has many solutions and to each solution A = Aj there 
corresponds a whole class of Einstein spaces containing a maximally symmetric space- 
time, being dS for Aj > 0, AdS for A^ < or Minkowski space (A^) for Aj = 0. For some 
Aj the maximally symmetric space may be stable. Each stable ground state (vacuum) 
defines a separate dynamical sector of the theory. Multiplicity of vacua for a L = f{R) 
gravity was first noticed in [TTj . 



We view (2) as an algebraic equation and assume that it has at most the countable 
number of solutions. Not every function f{R) admits a solution to (2). First note the 
degenerate case where any value of A is a solution (uncountable number of solutions); 
this occurs when (2) is viewed as a differential equation for /. Then f{R) = aR^ for any 
constant a 7^ [2S]- In the following we make some comments on this degenerate case. 
The eq. (2) has no solutions if its LHS defines a function of R, Rf\R) — 2f{R) = F{R), 
which nowhere vanishes. Treating this definition as a differential equation for / (for a 
given F) one finds that any Lagrangian which admits no CGS solutions is of the form 

f{R)=aR'±R'J^dR (3) 

with arbitrary F{R) > everywhere. Examples. 
1. Let F{R) = Ce^^ with 6, C ^ 0. Then 



f{R) =aR^ + C 



V^(l + 6i?) + -i?2 —AR 
2 ^ ' 2 J R 



the latter integral is non-elementary. 
2. For 



F{R) = Co + £ C2nR 



2n 

n=l 



*Any Lagrangian is determined up to a divergence of a vector field made up of the dynamical 
variables. If the gravitational Lagrangian is to be a scalar function of the Riemann tensor invariants 
alone and involve no derivatives of the curvature, the Lagrangian is determined up to a constant 
multiplicative factor. The factor must be fixed if any matter is minimally coupled to gravity in JF. 
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with Co and C2n > one gets 

m) = + aR' + c^R' In + ^ E -^R'""- 

^ ^ n=2 ^ 

Any gravity theory with a Lagrangian of the form (3) is unphysical and should be re- 
jected. Clearly there are infinitely many functions f{R) admitting solutions to (2) and 
thus possibly possessing a stable vacuum. A few examples of these Lagrangians. 

1. f{R) = R + aR^ + a-^R^, a > 0. There are three CGS solutions: M (Ai = 0), 
dS (A2 = +a) and AdS (A3 = -a). 

2. f{R) = -2A + R + aR'^ + bR^. 

Eq. (2) becomes a cubic equation A^ — A/6 + AA/b = having 3 distinct real roots if 
4A^ < (276)"^ and one real root Ai (and A2 = A3 complex) if 4A^ > (276)"^ In the 
limiting case 4A^ = (276)^^ there are 3 real roots with Ai = A2 and A3 7^ Ai. The 
solutions are independent of the coefficient a. In the case 6 = there is a unique CGS 
solution A = 4A, being dS or AdS, the same as in Einstein theory. 

3. For f{R) = -e""^, a > 0, there is a dS space with A = 2/a as a unique CGS solution. 

4. For f{R) = a\/R + bR'^^^, a, 6 > 0, there are two CGS solutions: Ai and 

A2 = +(i)^/-. 

The last example illustrates a general rule: if /(O) = then A = is always a solu- 
tion of (2), i.e., is a CGS solution even if /'(O) is divergent. In fact, if /'(O) ±00 
the term Rf'{R) may a priori either vanish at i? = or diverge either logarithmically 
or as inverse of a power law. If Rf\R) = for some n > then f = c — 
while for Rf{R) = (Ini?)", n > 0, one finds / = c + ;^(lni?)"+i; in both cases /(O) 
is divergent. For /(O) = and /'(O) — > ±00 the leading term in / near R — is R"', 
< a < 1, and then Rf'{R) — * 0. However if /'(O) and/or /"(O) is divergent the method 
for establishing whether A4 is stable does not work. 

We emphasize that in order to investigate the dynamics of a restricted NLG theory 
one needs exact solutions of eq. (2). We shall see that stabihty of a CGS solution is 
determined by the values of /'(A) and /"(A). In principle to check stability it is sufficient 
to find numerically an approximate solution A to eq. (2) and then approximate values 
of /'(A) and /"(A). Also the mass of the scalar component of gravity is determined 
by these two numbers. However an exact solution is necessary to calculate the scalar 
field potential both in Helmholtz-Jordan and Einstein frames; otherwise one gets only 
approximate equations of motion in these frames as is shown in the following example: 
L = f{R) = a sin a > 0. Introducing a dimensionless quantity x = j one finds 
that eq. (2) reads x cos x + 2 sin x = 0. cos a; = is not a solution and the equation 
may be written as a; -|- 2tga; = 0. The obvious root is a; = 0, but it corresponds to 
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R = X = oo and this solution must be rejected on physical grounds. In the interval 
— 7r/2 < x < 7r/2 where tangens is continuous the functions x and tgx are of the same 
sign and the equation has no solutions. In each interval {n — l/2)7r < x < {n + l/2)7r, 
n = ±1, ±2, . . ., the equation has exactly one solution which may be determined nu- 
merically. The scalar component of gravity is defined as p = df /dR and to determine 
the potential for p one needs to invert this relation to get R = r{p). In the present 
example p = —{a/RYcosa/R and though this relation is in principle invertible (since 
f"{R) 7^ and /" vanishes only at separate points where tga/R = 2R/a), it cannot be 
inverted analytically in any of the intervals. One sees that exact solvability of eq. (2) is 
often correlated to exact invertibility of the definition p = f'{R)- We conclude that the 
condition of exact analytic solvability of eq. (2) is of crucial importance and in practice 
imposes stringent restrictions on the Lagrangians excluding many simple combinations 
of elementary functions. A further constraint will be imposed in the next section. 

Finally we make two remarks on the field equations (1). 
Firstly, recall that for cosmologists the most attractive Lagrangians are those containing 
inverse powers of R rather than being polynomials in R. In consequence the coefficients 
of fourth order derivatives in (1) are rational functions and this implies that one should 
deal with great care with various terms in these equations in order to avoid multiplying 
or dividing by zercEl. For simplicity we demonstrate it on a toy model. Suppose that 
field equations read 

R^, + -^OR^u = 0. (4) 

Multiplying them by R^ one gets 

R'^R^, + UR^, = (5) 

and a class of solutions to these equations is given by R^^, = ip^y ^ where the tensor 
is traceless, R = ip = g^^ip^u = and satisfies Oip^^ = 0. However is not a solution 
to (4) since the LHS of these equations is then ip^^i, + 0/0. A class of solutions to (4) 
is of the form R^i, = (p^i, ^ and D^^;^ = —4>^4>\i.v with = g'^'^cp^i,; clearly these are 
also solutions to (5). Furthermore, any spacetime satisfying R^i, = is a solution to 
both (4) and (5). At first sight this is not since the second term in (4) becomes diver- 
gent. One may however give a precise meaning to this term by trying an Einstein space, 

^Wc stress that this is not trivial. In a frequently quoted paper [25] the trace of eqs. (1) for a 
Lagrangian R— 1/R was multiplied by R'^ giving rise to a scalar equation for R admitting i? = as a 
solution and thus Minkowski space; further considerations of the work were based on perturbations of 
this spacetime. Actually the field equations for this Lagrangian have only dS and AdS spaces as CGS 
solutions. This "curvature instability" found in [28] has been generalized to many other functions f{R) 
without checking if Minkowski space is a solution and is even regarded as advantage of the Palatini 
formalism over the purely metric gravity theories [29] . This error of introducing or omitting some 
classes of solutions by multiplying the field equations by a power of R may be traced back to Bicknell 
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Rflu = {^/4:)gfj.u, then D-R^i, = and the eqs. (4) reduce to Xg^^^, = so that i?^,^ = 
actually are solutions. In conclusion, by replacing the correct equations (4) by allegedly 
equivalent equations (5) one introduces a class of false solutions R^u = ip^y. 

Secondly, we comment on the cosmological constant |31,J. In metric NLG theories 
this notion has a rather limited sense. In general relativity A is both the constant 
appearing in the Einstein-Hilbert Lagrangian, A = — iL(0), and the curvature of the 
unique maximally symmetric ground state, A = i?/4. If /(O) 7^ is finite in an NLG 
theory one may define A as — 1/(0), however there is at least one COS solution with 
the curvature = A 7^ whose value is independent of the value /(O) (in the sense 
that the function F[R) = Rf'{R) — 2f{R) may be freely varied near R = provided 
-F(O) 7^ is preserved, then R = X remains the solution of (2)). Alternatively, A may 
be defined as A/4 for each vacuum (stable ground state), then A has different values in 
different sectors of the theory. However this cosmological constant is related solely to 
the vacuum and does not appear as a parameter in other solutions to the field equations 
(1). We therefore shall not use this notion. 



4 The field equations including matter 



We shall now express the field equations in the form appropriate for investigating stabil- 
ity of the CGS solutions. Detailed calculations based on the general formalism [221 133] 
are given in [34|. The scalar component of the gravitational doublet is defined in HJF 
as p = J^, this canonical momentum is dimensionless. The definition is inverted to give 
the curvature scalar i? as a function of p, R{g) = r{p), i.e.. 

The inverse function r{p) exists iff f"{R) 7^ 0. The pure gravity Helmholtz action 



Shj = J d'^x Lnig, p) 



with Lh = p[R{g) — r{p)] + f{r{p)) (see Paper I) gives rise to the field equations 

'1 



p D 



'-f{r{p)) + r{p) 



P 



and 



2 1 

°P - ^f{r{p)) + -pr{p) = 0. 



(6) 



(7) 



By taking trace of (6) and employing (7) one recovers the relation R{g) = r{p). The 
effective energy-momentum tensor for p contains a linear term signalling that the energy 
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density is indefinite and deceptively suggesting that all solutions, including the CGS 
ones, are unstable [B]. However 9^i, turns out unreliable in this respect and to study 
stability one makes the transformation from HJF to Einstein frame being a mere change 
of the dynamical variables. It consists of a conformal map of the metric, 

Qfj.i' * 9^ll' = P9^lv•l 

and a redefinition of the scalar, 

P = exp i^J^'^'^ 

or (f) = \J\\ Inp, with k being a dimensional constant to be specified lateil§. Under the 
transformation of the variables the action integrals in HJF and EF are equal, 

Shj = Se = j d^x ^f^LH{g,p{(t))), (8) 

what defines Lh- To get the total Lagrangian precisely as in general relativity one 
introduces an equivalent Lagrangian proportional to Lh, 

and sets (2k^c)~^ = c^/{16ttG) or /t^ = SnG/c^. Hence is a minimally coupled scalar 
field with a self-interaction potential, 

^-lr^<p,a<pp-v{pm. (10) 

The constant k, determines the dimension of 0, [0] = g^/^cm^^^s^^, while V acquires 
dimensionality of energy density. The field equations following from (9) are 

Gf,u{g) = K^Tf,y{(p, g) = /t^(0,^0,^ - ^g^,ug"^(p,a^,p - gf,uV{^)) (11) 

and 

Solutions for a self-interacting scalar field in general relativity were studied in many pa- 
pers, however they are not solutions to eqs. (10)-(12) since the potential (10) is in most 
cases different from the potentials appearing in those papers. For example, an exponen- 
tial potential Voexp(— a/t^) with constant a was investigated in a number of works (see 

®In Paper I, for simplicity we have put k = 1 in eqs. (8), (9) of that paper and the definition of (f>. 



^0 = -7^9 



1 

2^ 



rip) f{r{p)) 



P 
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e.g. [33); in terms of the scalar p it reads Vop~", but there are no simple Lagrangians 
L = f{R) generating this potential via eq. (10). Recall that as long as one considers 
pure gravity, i.e. there is no minimally coupled matter in JF, the original Lagrangian 
L = f{R) is determined up to an arbitrary constant factor A. Let f{R) = Af\R). Then 
p = f'{R) = Ap and the inverse relation is R{g) = f{p). On the other hand R{g) = r{p) 
so that f{p) = r{p) = r{p/A). This implies LH{g,p) = ALH{g,p), the conformal factor 
p generates in EF the metric g^i, = Ag^^ and LH{g,p) = A~^LH{g,p)- 
The conformal map should not alter the signature of the metric, thus one requires p > 0. 
In general f'{R) cannot be positive for all R and it is sufficient to require that the map 
preserve the signature at the CGS solutions, i.e., p(A) = /'(A) > for each solution of 
eq. (2). Then p > in some neighbourhood of i? = A. If p(A) < one should take the 
Lagrangian L = —f{R). It may occur for some f{R) having multiple solutions of (2) 
that p(Aj) > and p{Xj) < for i 7^ j, then one should appropriately choose the sign 
of L at each sector of the theory separately. We shall assume that this has been don^ 
and piXi) = fiXi) > 0. 

The transformation from HJF to Einstein frame exists in a neighbourhood of a CGS 
solution with i? = A iff /'(A) 7^ 0. If /'(A) = the EF does not exist and the method 
of checking stability of the CGS solution does not apply. From A/' (A) — 2/ (A) = it 
follows that /(A) = and assuming that / is analytic around i? = A it has a general 
form 

00 

fiR) = J2ar,iR-\r (13) 

n=2 

for any real A. Notice that the degenerate Lagrangian L = R^ belongs to this class. 
This class of singular Lagrangians needs a separate treatment (see section 6) and we 
assume that f{R) is not of the form (13). 

For Lagrangians which are different from (13) the potential V^(0) in EF is not a constant. 
To prove it one assumes that V = const and determines the corresponding f{R). From 
(10) one gets 

f(r) 

r{p) =Cp+ ^ (14) 

where C = 2k,^V and one differentiates this equation with respect to / employing 

dr /df\"^ 1 



One finds 



df \dr J p 

dr _ 1 _ dp 1 f dp dp f. 

df p df p p^df df p^- 



''One may try a simplificatioii by choosing L{R) = jtj^, then L'{X) = 1. Actually this choice does 
not simplify the expressions for derivatives of the potential V and we shall not apply it. 
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dp / d/ dr \ ^ 

Since — = — 7^ this yields 

a J \ar dp J 

f{r{p)) = Cp\ 

Inserting this value of / into (14) yields r{p) = 2Cp and substituting P = ^ from the 
latter relation back to / = Cp"^ one finally finds / = Using R{g) = r{p) one arrives 
at f{R) = f for any real C 7^ 0, i.e., the degenerate Lagrangian. In particular the 
potential cannot vanish identically. In fact, V = implies r{p) = f{r)/p. Differentiating 
this relation with respect to r under assumption that f'{r) 7^ and f"{r) 7^ (the con- 
dition for r{p) to exist) one arrives at ff'/p"^ = implying /" = 0. This contradiction 
shows that V ^ 0. 

For admissible Lagrangians the potential is variable and this feature will be used to 
establish stability. 

Finally we comment on stability of a CGS solution in the presence of some matter. 
In our opinion it is the stability of pure gravity (only the metric and the scalar) that is 
crucial for physical viability of the theory while exotic forms of matter violating DEC 
can make the ground state unstable even in general relativity as it found in the two 
examples mentioned in section 2. Yet recently there appeared claims (see e.g. [HS] ) 
that the very presence matter (perfect fluid stars) renders f{R) gravity unstable. We 
show now that this is not the case. The point is that the property of the variational 
matter energy-momentum tensor (the stress tensor for short) to satisfy DEC is preserved 
under a conformal map of the metric. If one assumes that Jordan frame is physical and 
minimally couples a given species of matter in this frame, then the stress tensor 
t^j,(\E', (7)satisfies DEC by assumption. The field equations (1) for the metric read therj^ 

f'{R)R^. - ^f{R)g,u - V,VJ\R) + g,.nf\R) = g). (15) 

The conformal map g^p = pg^v makes the matter Lagrangian explicitly dependent on 
the scalar gravity p and the stress tensor in EF for ^ alone cannot be unambiguously 
derived from it. It is therefore convenient to express the gravitational field equations in 
both the frames in terms of t^j^ which is already defined as the variational one in terms 
of the physical metric (i.e. in JF). The metric field equations in EF replacing (11) are 
then [M] 

G^M = K%M,~g) + -V(^> -9)- (16) 



^The choice of a constant coefficient in front of f{R) in the Lagrangian should be determined by a 
Newtonian Hmit of the theory. This can be unambiguously done in the case where Minkowski space 
is the ground state. If the CGS solution under consideration is dS or AdS, the Newtonian limit is not 
well defined and the coefficient is undetermined. This trouble does not affect the present argument. 
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Since p > in a vicinity of the CGS solution and DEC holds for both the stress tensors 
in EFIj, it also holds for the total stress tensor T^^ + This means that matter 

cannot destroy stabihty of the ground state if it is stable in pure gravity theory. We 
comment on the instability found in [36] in section 7. 



5 Positive energy theorem for anti— de Sitter space 

We emphasize that the applied here method of proving stability of dS, AdS or Jvi spaces 
is based on the assumption that the scalar component of gravity satisfies in EF the domi- 
nant energy condition, what is equivalent to > 0. The fact that it implies positivity 
of total ADM or AD energy is not used. Nevertheless we shall consider for the moment 
this energy. In [3l] we proved that if L = f{R) is analytic at = and its expansion is 
L = R + aR^ + . . . and the potential V{(j)) in EF is non-negative, the ADM energy of a 
spacetime which is asymptotically flat is the same in both Jordan and Einstein frames 
and is non-negative. Near Ai the potential behaves as V = \aR^ + 0{R^), whence 
> for a > 0. An analogous positive-energy theorem may be proved in restricted 
NLG theories for spacetimes which are asymptotically AdS space. The case of space- 
times which asymptotically converge to de Sitter space is more complicated because dS 
is not globally static and we disregard it. 

Let Qfj^iy be the metric of AdS space in the following coordinates: 

2 2 

= g.^dx^'dx" = -(1 + ^) dt^ + (1 + ^)~Mr2 + r\d9^ + sm^Od^^), (17) 

the cosmological constant is A = — a = const > and R = \ = 4A. Let g^^ be a 
solution of the field equations (1) in JF which asymptotically approaches AdS metric 
(17), Qfj^u = g^u + h^y. Clearly g^j,,, is a solution to Einstein field equations G^y{g) = 6^^ 
in HJF, then the Abbott-Deser approach [7] applies and the total energy of the fields 
g^u and p is given by their formula, which in the case of (17) reduces to 



EadIq] = -—^\im smeded(f)[-r^dihoQ + —dihn + ^{d2hu + ^^d3hi3) 
IGtcG °° J sm 6 

+3rhoo + ^r^'hi - ^(/i22 + + ^/ii2 ctg^], (18) 



here = (r, 6, (p) and the timelike Killing vector in the Abbott-Deser formula is chosen 
as C,'^ = 6q, then its normalization at r = is C.'^C,^ = —1. In general all the components 



^This frame invariance of DEC does not apply to the scalar gravity p since in any frame its stress 
tensor is defined as that equal in pure gravity (no matter) to the Einstein tensor of the given metric. 
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of h^i, are algebraically independent and the requirement that separately each term in 
the integrand of (18) gives rise to a finite integral (what amounts to requiring that each 
term be independent of r) provides the asymptotic behaviour of: 
hoo, /i22 and h^s are of order r~^, hu = 0(r~^) and = 0(r~^) = /113. 
A spacetime being asymptotically anti-de Sitter space is defined in [37j and according 
to this definition a solution approaches AdS slower than is required by finiteness of its 
energy. We assume that the six components of h^^, behave as shown above while the re- 
maining four components, which do not enter the energy integral, tend to AdS as in the 
definition in [37j, /loi = 0(r~^) and /102, ^03 and are 0(r). Under these assumptions 
the scalar R{g) for a solution with finite energy approaches R = 4A as i? ^ 4A + 0(r~^). 

In Einstein frame an analogous integral expression for EadIq] holds for the corre- 
sponding solution g^j^i, with h^j^^, replaced by h^^i, = p{R)h^y. For r ^ 00 the conformal 
factor is j9 = f\R) = f\R + 0{r-^)) = f'{AA) + 0{r-^) (assuming that /"(4A) ^ and 
finite), whence EAnid] = f'{4:A)EAD[g] is finite. This energy is positive according to the 
positive energy theorem in general relativity provided V{(j)) > 0. Since /'(4A) > by 
assumption, we get that in spite of the indefiniteness of the tensor On,^{g,p) in HJF the 
positive-energy theorem for restricted NLG theories holds: 

(i) if L = f{R) admits AdS space with ^ = 4A < as a solution, (ii) /'(4A) > and 
/"(4A) 7^ is finite, (iii) the potential V{(f)) in EF is non-negative and (iv) a solution 
g^j_y in JF or equivalently the pair {g^u,p) in HJF tends sufficiently quickly to AdS space 
for r ^ 00, then the total energy in JF is equal to the AD energy in HJF and positive 
and proportional to that in EF, 

EAD[g] = {f'mr'EAD[g]>0. 

Recall that the AD definition of conserved energy only makes sense in HJF (and EF) 
since we have no notion of total energy for fourth-order equations of motion. Total 
gravitational energy in Jordan frame is therefore defined as a quantity equal to that in 
HJF. 

6 Minimum of the potential and stability 

In order to establish whether the potential for the scalar gravity (p in Einstein frame 
is non-negative in a vicinity of a candidate ground state solution A^, dS or AdS, it is 
necessary to calculate the first and second derivative of V at this state. To this end one 
first determines the derivative ^ of the inverse function R = r{p) to the definition of 
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the scalar, P = It is equal to 



dR 
dp 



' dp^ 
dr . 











di?2 


R=r(p)_ 


r=r(p) 



(19) 



Applying this outcome to the potential in (10) one finds 



dV 



-f {r{p)) - r{p) 
P 



dp 2n'^p'^ 

and this expression should also be inserted into the field equation (12) for (p. 



(20) 



Consider a CGS solution in Jordan frame with Gap — — |A^a/3 and Rijj) = A where 
A is a solution to (2). In HJF the scalar p at this state is pq = p{\) = /'(A) > 0. For 
the function r{p) one has r{po) = r(/'(A)) = A. Under the conformal map from HJF to 
EF the metric g^i, of the CGS solution is mapped to g^i, = poQui, = f'{}^)g^u and the 



scalar is equal to 



3 1 



Mn/'(A). The Einstein tensor remains invariant under a 



constant rescaling of the metric, hence 
and this allows one to define a cosmological constant in Einstein frame as 



4/'(A) 



Thus A4, dS and AdS spaces in JF (and HJF) are respectively mapped onto M., dS and 
AdS spaces in EF satisfying G^uisi) = —^(jfiu and being the CCS solutions in Einstein 
frame. Physical excitations of the field (f) in EF should be counted from its ground value 



bo, i.e., are equal ip 



0, then p = /'(A) exp{J^Kip). The potential V at 



t>Q IS 



1 

2/^p2i 



^('/'o) = 7ct[A/'(A)-/(A)] 



and applying (2) it equals to 



The potential for the scalar excitation ip is then 

[/(^) ^ vi4>) - v(0o) = Vipi4>)) 



A 



(21) 
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and vanishes for vanishing excitation, U{0) = 0. The field equation (11) is now modified 
to (hereafter g^^^ denotes any dynamical metric in EF, not only the maximally symmetic 
CGS solutions) 



G^,{g) + Ag 



K 



The first derivative of U with respect to ip (or p) vanishes when ^f{r{p)) — r{p) 
and this equation viewed as an equation for r coincides with eq. (2). Hence ^ 



(22) 


only at the CGS solutions with r{pi) = Aj = r(/'(Aj)), i = 1, . . . ,n. In other terms the 
equation A/' (A) = 2/ (A) determines all stationary points of U. At each of these points 
the potential Uilip) = V{(f)) — Aj/k^ vanishes provided Aj = Aj(4/'(Aj))~"^. On the other 
hand U (ip) (for a fixed value of A) may also vanish at some points different from the 
solutions Aj but these are not its stationary points; if such points exist the dominant 
energy condition is broken and a kind of (nonlinear) instability may develop. 

The second derivative of the potential, determining its behaviour at a stationary 
point is, from (20) and (19), 



dV^2 



1 

3^ 



-4/(r) + pr 



p^ 



f"{r)_ 

At the stationary point R = r = \, tp = and p = /'(A), then 

d^U -A 1 



+ 



d^|J^ ^=0 3/'(A) 3/"(A) 



(23) 



(24) 



For regular Lagrangians we are considering in this work one has /"(A) ^ finite. The 
potential U {ip) is non-negative if it attains minimum at ^ = 0, i.e. 



d^U 



di)^ 



V'=o 



> 0. 



Stability conditions were also derived by other authors applying different methods. Our 
condition is equivalent to that found in [23j which after using A/'(A) = 2/(A) reads 



3/'(A) 



d^U 



d^/^2 



V'=o 



> 0. 



(25) 



The linear perturbation method applied in [22] implies that stability occurs whenever 
the weak inequality in (25) holds. Cognola et al. [38] employ a minisuperspace approach 
to the stability problem (perturbations are spatially homogeneous) and get a stability 
condition of de Sitter space which is equivalent!^ to ours; also computing one-loop 
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The formula for the condition seems to be misprinted since it disagrees with their earUer work. 
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quantization corrections to L = f{R) they find this condition for dS background 
Yet Song et ah [101 UT] define stabihty of spatially fiat R-W spacetime in a nonstandard 
way: a gravity theory is stable if it approaches general relativity at high curvatures (for 
small R the theory should diverge from GR by definition); this cosmological criterion 
does not deal with a ground state solution. 

The derivation of (24) holds both for A = and A 7^ 0. The case A = is simpler to 
study. In this case /(O) = and assuming analyticity around R = one has 

00 

f{R) = R + aR"^ + J2anR'', (26) 

n=3 

then /(O) = 0, /'(O) = 1, /"(O) = 2a =^ and f/"(0) = ^. For a > the potential 
U > and the scalar field satisfies DEC. For spacetimes which are asymptotically fiat 
it is known [33] that EADM[g] = -E'^DAf[^, V'] ^ and the total energy vanishes only in 
Minkowski spacetime, g^u = Vtiy = 9tiy ^i^d ip = 

In de Sitter space (A > 0) there are following cases: 

— for /"(A) < the potential attains maximum at ■?/' = and the space is unstable; 

— for /"(A) > and /'(A) > A/"(A) there is minimum of U and U{'il)) > 0, hence the 
space is stable; 

— for /"(A) > and /'(A) < A/" (A) one finds U"{0) < and instabilities develop. 
For anti-de Sitter space the situation is reversed: 

— for /"(A) < and /'(A) > A/"(A) the negative potential attains maximum and the 
space is unstable; 

— for /"(A) < and /'(A) < A/" (A) the potential is at minimum and AdS is stable; 

— for /"(A) > the minimum of U shows stability of the space. 

Finally we return to the problem of singular Lagrangians (13) for which /(A) = = 
/'(A); for them the derivative (24) is divergent and the method of deriving it does not 
work. One may instead apply the gauge invariant perturbation method for de Sitter 
space directly in Jordan frame which gives rise [23] to the inequality (25). Let the lowest 
nonvanishing coefficient in the series (13) be a^. If > 2 then also /"(A) = and the 
expression (25) becomes indeterminate. In order to give it a definite value we define a 
function 

J(i?)^-i? + -^ (27) 

and define J(A) as its limit for _R — > A. Let i? = A + e, |e| <^ 1, then f'{R) = 
kaue^-^ + O(e^), f"{R) = k{k - l)ake''-^ + Oie''~^) and 

J(A + e) = -(A + e) + ^ + 0(e2). 
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Hence the stability criterion is J{X) = —A > 0. Recall that the method works only in dS 
space, A > 0, therefore the conclusion is that for all NLG theories having Lagrangians of 
the form (13) with A > 0, de Sitter space (as a CGS solutioE0) is unstable. None of the 
methods can be applied to these Lagrangians in the case A < 0. It might be argued that 
by continuity the criterion J(A) > should also work for A < 0, then all AdS spaces 
would be stable in these theories. However this argument is of little reliability. 

In Paper I an astonishing theorem was mentioned to the effect that an anti-de Sitter 
space may be stable in spite of the fact that the scalar has a tachyonic mass (i.e., 
the potential U{ijj) < and attains maximum at this space) [12] ■ In fact, if small 
fluctuations of the scalar gravity vanish sufficiently fast at spatial infinity of AdS space 
(i.e., for r ^ cxD in the metric (17)), the kinetic energy of the field dominates over its 
negative potential energy and the total energy of the scalar. 



EW = - J d^xv^T^., 



where T''^ is given in (22), is finite and positive, < E{iIj) < oo. This occurs if > |A 
at tp = 0. Since the energy of gravitational perturbations of AdS space is positive [7], 
the total energy of metric and scalar field fiuctuations is positive and Breitenlohner 
and Freedman conclude [12] that AdS space is stable against these (small) fiuctuations. 
Applying the definition of A in Einstein frame arising in NLG theories, the condition of 
stability of AdS in the case of maximum of the potential reads 

d^U 3A 
° > dV^ ^=0 > 16/'(A)' ^ ^ 

It should be stressed, however, that in this case the DEC is violated (only the total energy 
of the scalar is positive). From the viewpoint of a rigorous mathematical approach to 
the stability problem the condition (28) is rather unreliable |20j . 



7 Examples: specific Lagrangians 

We now apply the stability criteria of the previous section to a number of Lagrangians, 
some of which were already discussed in the literature. We assume that the Lagrangians 
depend on one dimensional constant fi and some dimensionless constants, fi is positive 
and has dimension (length)"^ so that R/fi"^ is a pure number. 



4n+4 

1- L = R + ^^, n = 0,l,.... (29) 



^^Besides R — X there are in general other solutions to Rf'{R) — 2f{R), e.g. for f{R) = a(i? — A)^ 
the other solution is R — — 2A. 
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This Lagrangian belongs to the class which admits no CGS solutions since it is given 
by eq. (3) for a = and F{R) = -R - (2n + 3)/x^'^+^i?-(2"+i). Clearly it should be 
rejected. Yet according to Sawicki and Hu [H] the theory for n = converges to general 
relativity for large R and in this sense is admissible. 



2. L = R+i^, n = l,2,.... (30) 

There is only one CGS solution with A = A_ = -{2n + 2)^ fx'^ < 0, /'(A) = {2n + 
l){n + 1)^^ > and the scalar is 

p{r) = l-2niyj . (31) 

We consider spacetimes with i? = r in vicinity of i? = A_, so that — oo < r < and 
1 < p < +00. The inverse function and the potential are respectively 

1 



with A_ < 0. The potential is always non- negative and U < 



A_ 



It attains minimum 



at ip = showing that AdS space is a stable ground state solution for this theory. 



3. L = R + for - 1 < n < real, (33) 

for non-integer n the function i?" is replaced by |-R|". The background evolution of 
the R-W spacetime is cosmologically acceptable [2] and solutions of the linear per- 
turbation equations for this Lagrangian are not incompatible with the observational 
data |l3]. However the equation for a ground solution gives rise to the contradiction 
|^|n+i _ _|_ 2)/i2»i+2 iniplying that the theory in untenable. 



,An+2 

4. L = R-^, n = l,2,.... (34) 

Here A and r{p) have the same moduli and opposite sign to those in the case 2: dS 
space is the unique CGS solution for A = A+ = — A_ > and r{p) > 0. Accordingly, 
A = A+ = — A_. Now we take r around r = A+, and again 1 < p < oo. The po- 
tential is U = f/+ = — f/-, hence it is contained in the interval — ^ < t/+ < 0. This 
indicates that U has maximum at ip = and this fact is confirmed by a direct computa- 
tion. In conclusion, de Sitter space is unstable and this theory is discarded as unphysical. 
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4n+4 

5. L = R-^^, n = 0,l,.... (35) 

This Lagrangian has been most frequently studied in apphcations to the accelerating 
universe, usually for n = 0. Most expressions here are akin to respective ones in the 
case 4. The field p is always greater than 1 and there are two CGS solutions for X± = 
±{2n + 3)2^/i^, hence p(A±) = p(A_) = The two CGS solutions define two 

different sectors of the theory which should be separately studied. 
A. De Sitter space sector. 

A = A+ > and the sector comprises all positive values of r. The inverse function is 

r{p) = r+(p) = (^^-^l (36) 



giving rise to the potential 



2n+l 



C(p(,») = f/, = !i±I(Z^)-«,^-^ (37) 

which is always non-positive and attains maximum at dS space. This space is then 
unstable (for = it was found in [HI |23l |39]) and this sector of the theory must be 
rejected (on other grounds this conclusion was derived in |45j). 

Seifert [36j finds that gravity theory (35) is highly unstable in the presence of matter: 
a static spherically symmetric solution becomes unstable to linear spherically symmet- 
ric perturbations if perfect fluid matter forms a quasi-Newtonian polytropic star. This 
result is derived applying an intricate variational method and requires very long com- 
putations. We note that (besides the fact that the Newtonian limit is not well defined 
there) the author assumes that R is approximately equal to the stellar matter density. 
This means that he deals with spherically symmetric perturbations of dS space. Since 
this space is unstable in pure gravity (35) it would be rather surprising if a small amount 
of matter could stabilize it. 



B. Anti-de Sitter space sector. 
Its existence (for n = 0) was first noticed in ^3], then in |l6], but its properties were 
never analyzed in detail, probably due to the fact that a negative A does not fit the 
observed accelerated expansion. A = A_ < and accordingly — oo < r < 0, hence 
r{p) = r_{p) = —r^{p) and U = = — 1/+ with A_ = — A_|_. This potential is non- 
negative and has minimum at ip = 0. This sector has a stable ground state solutioij^ 
and in this sense it forms a viable gravity theory. The scalar gravity has mass being 

"'^^In [Uj it is claimed that Lagrangians given in cases 4 and 5 (for both n even and odd) always 
develop instabilities while Lagrangians in cases 1 and 2 always describe a stable theory. 
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a function of n, for n = it is = -^j^/i^ while for n — oo it tends to 
Disregarding the incompatibihty of this theory with the cosmic acceleration, one may 
make a rough estimate of /i. Since A is of order — for all n > and the observational 
limit is |A| < lO^^^m"^ one gets an upper limit /i < 10~^^m~^ or fihc < 10~^^eV, very 
small indeed. 




6. L = In — , g > 0. (38) 



One assumes -R > 0. A unique possible ground state is dS space with A = 6*^/1^ , then 
p{X) = 2exp(glng) and /"(A) = /i~^(2g — 1) exp[(g — 1) Ing — g]. The Lagrangian must 
be regular, i.e. /"(A) ^ implying q ^ 1/2. From (24) 

U"(0) = ^ e'ja-ing) 

^ ' 6(2g-l) 

hence for < g < 1/2 dS space is unstable while for g > 1/2 de Sitter space is a stable 
ground state. The function p{r) may be inverted and then the potential can be explicitly 
calculated only for g = 1 or 2. According to [2] this theory is cosmologically acceptable 
for any g > though the matter era begins too early and its duration is too long. 



7. L = aR ^— ^, a>0. (39) 

sinh4 ~ 

This Lagrangian appeared in the metric-affine approach to gravity [48J. The equation 
Rf'{R) — 2f{R) = cannot be analytically solved even in the case a = (it can only 
be shown that the roots do not lie close to i? = 0) and for practical reasons this theory 
must be rejected. 



8. L = ^i'(\n^ + l]+^R\ a>0. (40) 



One may start from a more general Lagrangian 

L = -fR + b\n(cR) + ^R^, a',6,c>0, 

but then eq. (2) for A cannot be solved analytically. We therefore set 7 = and 
multiply L by fi'^/b and define a = y/i^; finally we choose such value of c as to get a 
simple expression for A. A unique solution to (2) is then A = /x^ and 

/^^ 2a 

p=— + —r, 41 
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r > 0. To invert this function we first notice that p{r) — > oo for both r — > and r — oo 
and has minimum at Tq — fj,^/y/2a equal to p(ro) = 2y/2a. Hence p{r) may be inverted 
either in the interval < r < ro or r > tq. To choose the correct interval one must 
establish whether A = /i^ belongs to the ascending or descending branch of p{r) and this 
depends on the value of a. We assume a > 1/2, then /i^ > tq and dS space lies on the 
ascending branch of p (for a < 1/2 a similar procedure can be performed). Solving (41) 
one chooses the larger root (both roots are positive) , 

rip) = ^{P + yjp^ - 8«) (42) 



since r — > oo corresponds to p — > cx). The potential is 



u ' 



P{p) - ^(lnP-ln(4a)) 



- ^, (43) 



where P = p + ^/p^-8a and A = This implies /'(A) = 2a + l > A/"(A) = 2a-l 

and the potential has minimum at -0 = 0. This theory has dS space as a stable ground 
state and is viable. 

The case a = 1/2 is singular since f"{n^) = and p{r) cannot be inverted around 
r = jji^ while /'{n"^) — 2. Formally the conformal map to EF exists at this point but the 
potential U cannot be defined there. None of the methods to check the stability does 
work there and it is reasonable to disregard this case. 

9. The limiting case a = of the Lagrangian (40) requires a separate treatment. 
Again \ = jj?, p = and r{p) = \j? jp > 0. /'(A) = 1 and /"(A) = — ^ give rise to 
U"{0) = — |a*^- De Sitter space is unstable making the theory untenable. 
The additive constant appearing in this Lagrangian (as well as in the case 8) is inessential 
in the sense that it only affects the absolute value of A (but not its sign) and has no 
influence on stability properties of dS space. In fact, for a Lagrangian 




a real dimensionless, one gets again p — and the value of A is shifted to A = 

//^exp(| — a), hence it is still dS space. Then U"(0) — — ^ [1 + exp(l — 2a)] imply- 
ing instability of the space for any a. This case is, however, exceptional: we will see 
below that in general not only /j, but also dimensionless parameters in L determine sta- 
bility of CGS solutions. 



10. L^i^'i-] (44) 
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for a rational (negative and positive) has also attracted some attention [23, SHI EH] since 
it is a scale-invariant theory. For non-integer a one takes |-R|". If a < the equation 
Rf'{R) = 2f{R) is solved only hj R = ±00 and we reject this case. For a = 2 one gets 
the degenerate Lagrangian R^ which we disregard. For a > 2 integer this is a singular 
Lagrangian (13) discussed in section 4 having A = and the criterion —A > yields 
that Minkowski space is stable for these theories. Putting aside the obvious case a = 1 
one considers a > non-integer. /(O) = always. For < a < 1 both /'(O) and 
/"(O) are infinite, for 1 < a < 2 there is /'(O) = and /"(O) = 00 and for a > 2 both 
/'(O) = /"(O) = 0. Once again one may apply the function J{R) defined in (27) and 
it is equal J = -R so that J(0) = and for all three cases the criterion J(0) > is 
satisfied. One may thus claim that for all a > Minkowski space is the unique stable 
ground state, nevertheless it is difficult to avoid impression that for a 7^ 1 the theory is 
bizarre and rather unphysical (and furthermore in conflict with the astronomical obser- 
vations, as mentioned in Paper I). 



a real [23 [211 [23] • There are two CGS solutions with A-t = ±v^/i^, which are the 
same as for the case a = (Lagrangian (35) for n = 0) since the R^ term does not 
contribute to A. The attempt to flnd the inverse function r(p) leads to a cubic equation 
and solving it would be impractical. We therefore quit from computing the explicit form 
of the potential (an implicit form of V is given in [2S]) and restrict ourselves to studying 
its extrema. 

A. De Sitter sector for A = A+. 

The condition p(A+) > requires a > —^rjn- This condition does not determine the sign 



and from (24) one flnds: 

— for — < a < dS space is unstable and 

— for a > dS space is stable. Yet cosmologically the theory in this case is inaccept- 
able since there is no standard matter era preceding the acceleration era 0. 

We omit the singular case a = where /"(A+) = 0. 
B. Anti-de Sitter sector with A = A_. 

Now the condition p(A_) > requires a < From (24) one gets that for — < a < 
the potential has minimum at ip = and AdS space is stable. Yet for a < —-^^ the 
potential has maximum. This, however, does not automatically imply the instability 
since one should furthermore apply the criterion (28) of positivity of scalar fleld energy. 
It follows from it that AdS space is 




(45) 



of 
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— stable for — < a < with respect to scalar field perturbations with positive 
energy, 

— unstable for a < — 

In the range of values of a for which the theory is stable in the standard sense (the 
potential has minimum) the mass of the scalar gravity excitations above dS space is 

ml = 3\/3[2(3V3a + 2)(3V3a - 1)]~V^ 

while in the case of AdS ground state it is 

m'i = 3V3[2(3V3a + 1)(2 - 3V3a)]-^fi^. 

The particle masses tend to infinity when a approaches the finite limits of the admissible 
range. m+ monotonically decreases and becomes very small for large values of a while in 
the AdS sector the scalar particle mass attains minimum = -^yU^ at a = (6v^)~^- 
In the interval < a < the theory has two viable sectors: one with dS space 
ground state for A+ = 3-\/3[8(3-\/3a + 2)]~^/i^ and the other having AdS as a ground 
state with A_ = — 3v^[8(2 — 3^/3a)]~^ jj? . Classically these are two different physical 
theories, each with a unique ground state. One cannot claim that this is one theory 
having two different and distant (in the space of solutions) local minima of energy. En- 
ergetically these two states are incomparable, each of them has vanishing energy (defined 
with respect to itself) and assuming that one of these minima is lower than the other 
is meaningless [5]. One may only compare the masses of the scalar gravity in the two 

theories. The mass ratio decreases monotonically from infinity for a approaching 

(3\/3)-i to zero for a tending to 2(3-\/3) ^ If one believes that this Lagrangian describes 
the physical reality a difficult problem arises: how does the nature choose which of the 
two theories with the same Lagrangian is to be realized? In our opinion the nature 
avoids this problem merely by avoiding this Lagrangian (and other ones with the same 
feature). 

This Lagrangian illustrates a general rule: all the parameters appearing in a Lagrangian 
do contribute to determination of stable sectors (i.e., physically distinct theories) corre- 
sponding to it. 



12. L = Rexp(^^^, e = ±l. (46) 

To each value of 6 there is one CGS solution with A = — ^//^ and p{X) = 2/e > 0, 
then U"{0) = —^6efi^. For 6 = —1 the potential attains minimum and de Sitter space is 
classically stable. For 6 = +1 one applies the stability criterion (28) for AdS space and 
one gets that this solution is unstable according to this condition too. The Lagrangian 
(46) for ^ = — 1 is also cosmo logically preferred since it is asymptotically equivalent to 
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the ACDM model [2\. Unfortunately the function p{r) cannot be inverted analytically 
and the explicit form of the potential is unavailable. 



13. Finally we consider a class of "toy models" ' possessing infinite number of ground 
states. For convenience we introduce a dimensionless variable x = R/ //^ and assume 

L = fiR) = f^'Fix) = /x^e^^^-) (47) 

where 

I{x) = f (48) 

^ ^ J x + h{x) ^ ' 

and h[x) is a continuous periodic function taking both positive and negative values, 
Ml <h{x)< M2 with Ml < and M2 > 0. The scalar field is 

V = ~r~ = 

ax X + h[x) 

and is positive if x + h{x) > 0. For an arbitrary h one cannot find r{p) and the potential; 
here it is sufficient to determine CGS solutions and U" at these states. The eq. (2) takes 
now the form x^ = 2F and since F > it is equivalent to 

ax 

On the other hand from the definitions (47) and (48) it follows that 

-^^ = x + hix), (50) 

dx 

hence those x which are solutions of (49) must also be solutions to h{x) = 0. Since 
Ml < h{x) < M2 there is at least one root of h{x) = and for a continuous periodic 
function there is infinite number of zeros, h{xn) = 0, n = 0, 1, . . . and A„ = Note 
that Xn 7^ since A = implies /(O) = /x^ exp(2J(0)) = while J(0) is finite by its 
definition. The function x + h{x) tends to ±00 for x —>■ ±00, hence there is a point 
X = y such that y + h{y) = and y 7^ 0. To ensure that x + h{x) > for x > y one 
requires x + h{x) be monotonic, i.e. 1 + h'{x) > 0. Then I{x) is defined (and positive) 
for all X > y. Denoting /„ = I{xn) one finds that U" at a point R = A„ is 

dV^ A„ 6 " l-h'iXnY ^ ' 

The condition h'{x) > — 1 does not determine the sign of the fraction and to this aim 
one must specify h. Here we choose as an example h{x) = ^{sinx — cos a;). Clearly 
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h'{x) = ^{sinx + cosx) > — 1 and the unique solution of a; + ^(sina; — cosx) = is 

y — 0, 3183 The zeros Xn > y oi h are solutions to tgx = 1 (cosx 7^ 0) and these 

are Xn — j+ nn, n = 0, 1, At these points h'{xn) — (— l)"^ and for n odd there is 

d'^ A2n+1 

therefore the infinite sequence of dS spaces with curvatures A2n+i = A*^3;2n+i defines 
unphysical (unstable) sectors of the theory. Yet the other sequence for n even consists 
of dS spaces having curvatures X^n = (2n + |)7r/x^ which are stable for this Lagrangian. 
The scalar particles corresponding to these sectors have masses 

^2 ^ (2n+|)^ 
6(V2-1) 



8 Conclusions 

In this paper wc have investigated stability of ground state solutions in L = f{R) grav- 
ity theories being cither Minkowski, de Sitter or anti-de Sitter spaces. Stability may be 
studied in any frame and Einstein frame is particularly suitable to this aim since one 
may apply there the methods developed in general relativity. We have given an explicit, 
effective and simple method of checking stability of these spaces based on the dominant 
energy condition applied to the scalar component of the gravitational doublet. After 
applying the method to thirteen specific Lagrangians (their ground states arc de Sitter 
and/or anti-dc Sitter spaces) corresponding to 20 different cases (depending on values 
of parameters in L) it was found that, as it was a priori expected, half of them give 
rise to viable theories (9 viable versus 11 untenable ones). And a generic feature is the 
existence of multiple vacua (stable ground states), each generating a separate physical 
sector or rather a separate gravity theory, all having the same Lagrangian. Hence it is 
expected that there is an infinity of viable gravity theories. What to do with such a 
wealth of theories (all differing from each other only by the form of the potential for the 
scalar gravity field)? 

We stress that it is incorrect merely to search for a theory which easily and immedi- 
ately accounts for the big problem of cosmology — the acceleration of the universe. After 
all general relativity was not formulated to solve some urgent problems in celestial me- 
chanics (the perihelion shift of Mercury) or in cosmology (non-existence of Newtonian 
cosmology) and for many years its confirmation was quite marginal. At the time of its 
advent its advantage was that it was physically much deeper and more general than 
Newton's gravity. And the same should be expected about a modified gravity which 
may ultimately replace Einstein's theory. Its physical content will be more relevant than 
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immediate observational confirmation. 

Before a deep creative pliysical idea will appear we need further viability criteria to 
maximally reduce the set of viable gravity theories. Undoubtedly one of the most im- 
portant ones will be the condition that a tenable theory must be in agreement with the 
Newtonian and post-Newtonian approximations to gravity — as soon as these approx- 
imations will be rigorously defined in de Sitter and anti-de Sitter backgrounds. One 
should, however, expect that the selected set will still be large and possibly infinite. 
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